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Although there is mounting observational evidence that the cosmic expansion is undergoing a 
late-time acceleration, the physical mechanism behind such a phenomenon is yet unknown. In 
this paper, we investigate a holographic dark energy (HDE) model with interaction between the 
components of the dark sector in the light of current cosmological observations. We use both the 
new gold sample of 182 type la supernovae (SNe la) and the 192 SNe la ESSENCE data, the baryon 
acoustic oscillation measurement from the Sloan Digital Sky Survey and the shift parameter from 
the three-year Wilkinson Microwave Anisotropy Probe data. In agreement with previous results, 
we show that these observations suggest a very weak coupling or even a noninteracting HDE. The 
phantom crossing behavior in the context of these scenarios is also briefly discussed. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The current idea of a negative-pressure dominated uni- 
verse seems to be inevitable in light of the impressive 
convergence of the recent observational results (see, e.g., 
0) S H & Hi)- This in turn has led cosmologists to hy- 
pothesize on the possible existence of an exotic dark com- 
ponent that not only would explain these experimental 
data but also would reconcile them with the inflationary 
flatness prediction (SlTotai = !)■ This extra component, 
or rather, its gravitational effects is thought of as the first 
observational piece of evidence for new physics beyond 
the domain of the standard model of particle physics and 
constitutes a link between cosmological observations and 
a fundamental theory of nature (for some dark energy 
models, see @. For recent reviews, see also Q)- 

On the other hand, based on the effective local quan- 
tum field theories, the authors of Ref. [8( proposed a 
relationship between the ultraviolet (UV) and the in- 
frared (IR) cutoffs due to the limit set by the formation 
of a black hole (BH). The UV-IR relationship in turn 
gives an upper bound on the zero point energy density 
Pa < AlpL~ 2 , which means that the maximum entropy is 

of the order of Sg^ . This zero point energy density has 
the same order of magnitude as the dark energy density 
@ , and is widely referred to as the holographic dark en- 
ergy (HDE) [Si (see also [U). However, the HDE model 
based on the Hubble scale as the IR cutoff seems not to 
be able of leading to an accelerating universe [3] . A so- 
lution to this matter was subsequently given in Ref. [13] 
that discussed the possibilities of the particle and the 
event horizons as the IR cutoff, and found that only the 
event horizon identified as the IR cutoff gave a viable 



HDE model [13]. The HDE model based on the event 
horizon as the IR cutoff was found to be consistent with 
the observational data [12f . 

A subsequent development concerning the idea of a 
holographic dark energy is the possibility of considering 
interaction between this latter component and the dark 
matter in the context of a holographic dark energy model 
with the event horizon as the IR cutoff. As an interesting 
result, it was shown that the interacting HDE model re- 
alized the phantom crossing behavior [l3|], which is also 
obtained in the context of non-minimally coupled scalar 
fields (see, e.g, [lj] and references therein). Other recent 
discussions on interacting HDE models can be found in 

©Gap. 

In this paper, we test the viability of the interacting 
HDE model discussed in Ref. [13] by using the new 182 
gold supernovae la (SNe la) data 0, the 192 ESSENCE 
SNe la data [3], the baryon acoustic oscillation (BAO) 
measurement from the Sloan Digital Sky Survey (SDSS) 
[H, and the shift parameter determined from the three- 
year Wilkinson Microwave Anisotropy Probe (WMAP3) 
data [j|. We organized this paper as follows. In Sec. 
II, we review the basic equations of the interacting HDE 
model considered in this paper. The choices of the IR 
cutoff are also discussed in this Sec. II. We fit the model 
based on the event horizon as the IR cutoff to the ob- 
servational data above mentioned in Sec III. We end this 
paper by summarizing our main conclusions in Sec. IV. 



II. INTERACTING HDE MODEL 

We consider a spatially flat Friedmann-Robertson- 
Walker universe with dark matter, HDE and radiation. 
Due to the interaction between the two dark components, 
the balance equations between them can be written as 
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FIG. 1: The results of our joint analysis involving the ESSENCE (192 SNe la) plus BAO plus CMB shift parameter, a) 
Confidence contours (la, 2a and 3a) in the b 2 — c parametric space. As discussed in the text (see also Table I), at 68.3% c.l. 
we find c = 0.851q'q2, and b 2 — 0.002+g 002 • b) Similar results for the c — Q m o plane. 



where the HDE density is 

PD = 3c 2 M 2 L~ 2 



(3) 



In the above equations, L is the IR cutoff, M p = 1/y/SwG 
is the reduced Planck mass, lod is the equation of state 
of the HDE, T = 9b 2 M 2 H 3 is a particular interacting 
term with the coupling constant b 2 , and the subscript 
means the current value of the variable. The HDE, dark 
matter and radiation density parameters are defined, re- 
spectively, as n D = p D /(3H 2 M 2 ), n m = p m /(3H 2 M 2 ), 
and f2 7 = Py/(3H 2 M 2 ). Note that, if we choose the 
Hubble scale as the IR cutoff, i.e., L = 1/H, then we 
find that fi m /Slc = (1 — c 2 )/c 2 , which means that the 
HDE always follows the dark matter. Even though the 
HDE equation of state wd can be less than —1/3 with 
the help of the interaction [l|| , this model cannot explain 
the transition from deceleration to acceleration. 

As suggested by Li [l(| , one can choose the future event 
horizon or particle horizon as the IR cutoff. For the fu- 
ture event horizon, 



L(t) = a(t) 



= a 



da' 



it o(f) J a H'a 
whereas for the particle horizon, 

, s t* dt ' r da ' 



o(f) 



(4) 



(5) 



Substituting Eqs. ([4]) and (J5J) into Eq. Q and taking 
the derivative with respect to x = In a, we obtain 



n L , = — = -6M 2 H 2 n D ± ^// 2 <>; ; 2 . 



, _ dp D 
dx 



(6) 



where the upper (lower) sign is for the event (particle) 
horizon. Since po = dp^/dt = p' D H, Eq. ^ can be 



written as 



p' D + 3(1 + lu d ) Pd = -9M 2 b 2 H 2 



(7) 



Combining Eqs. ([6]) and ([7]), we obtain the equation of 
state of this interacting holographic dark energy, i.e., 



1 2 
3^3 



(8) 



When the interaction is absent, b 2 = 0, it is clear from 
the above expression that we cannot choose the particle 
horizon as the IR cutoff. In 17|, it was argued that the 
effective equation of state of the HDE in the interacting 
case should be 



eff 

= UJ D 



(9) 



Based on this effective equation of state, it was concluded 
that there was no phantom crossover even for an inter- 
acting HDE model. In fact, by combining the Friedmann 
equation with Eqs. HJ and (|2|), we obtain the accelera- 
tion equation 



H = -4irG(p + p). 



(10) 



For a flat universe, the physical consequence of the phan- 
tom dark energy is a super-acceleration when the dark 
energy dominates. Note that it is ojd, not wff , that ap- 
pears in the acceleration equation (fT0|) . Therefore, the 
effective equation of state seems not to show the true 
physical meaning of the equation of state of the HDE, 
and ud should be used instead. 

Substituting Eq. (jHJ) into Eq. ((7]) and applying the 
definition of f2u, we have 
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FIG. 2: a) The same as in Figure 1 when the ESSENCE (192 SNe la) data are replaced by the new 182 Gold sample. The 
contours in the b 2 — c plane also correspond to la, 2a and 3<r. b) The same as in the previous Panel for the c — Q m o parametric 
plane. 



On the other hand, substituting H = H'H and pd = 
udPd into Eq. (|10|) . we obtain 



W_ 
JT 



2 c 



3/2 



1 



-fL 



(12) 
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Now, combining Eqs. (jlip and (| 1 2[) . we find the differen 
tial equation for f^, i.e., 



f2' n _ 2-v/f! / 



_£ = i - n D ± 



(1 - O^,) - 36 2 + IX, 



(13) 



a result that is consistent with Eq. (5) of Ref. [18| when 
the radiation term f2 7 is neglected. 

If we choose the particle horizon as the IR cutoff, the 
current acceleration requires that lujjq < —1/3 — (£l m o + 
2f2 7 o) /3fi do- From Eq. ||SJ), we also obtain a lower bound 
on b 2 , 



b 2 > 



n, 



9 fi 3 / 2 
z "do 



(14) 



The past deceleration and the transition from decelera- 
tion to acceleration requires that fl' D0 > 0, so Eq. JT 
gives the upper bound on b 2 , 



b 2 < 



n 



i»0 



1 - 2 



DO 



1 



-,{)■ 



(15) 



By comparing Eqs. (fl"4"l) and (fT5|) . we see that the upper 
bound is lower than the lower bound, so that the inequal- 
ities are not satisfied. The model based on the particle 
horizon as the IR cutoff is not, therefore, a viable dark en- 
ergy model. In what follows, we consider only the HDE 
based on the event horizon. As discussed in [1^, the 
interaction T cannot be too strong and the parameters 
b 2 and c are not totally free; they need to satisfy some 
constraints. Following [l8|, we take < b 2 < 0.2 and 
VTTd' <c< 1.255. 



III. OBSERVATIONAL CONSTRAINTS 
A. The data 

There are three parameters ft m o, c and b 2 in the in- 
teracting HDE model since Q jjq — 1 — f^mo — and 
f2 7 o ~ 1CP 5 . In order to place limits on them and test 
the viability of the model, we apply both the 182 gold 
SNe la and the 192 ESSENCE SNe la data [| to fit 
these parameters by minimizing 



E 



(16) 



where the extinction-corrected distance modulus fi(z) = 
51og 10 (di(2)/Mpc)+25, o~i is the total uncertainty in the 
[lobs observations, and the luminosity distance is given by 



(1 + *) 



dz 



H{z>) 
z) 2 



l F(z)V^W 



(1 + *) 



(17) 



where z = ao/a — 1. In all the subsequent analyses, we 
have marginalized the Hubble parameter Hq- 

In addition to the SNe la data, we also use the BAO 
measurement from the SDSS data [4, 5] 
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(18) 



and the CMB shift parameter measured from WMAP3 



4 











c=0.85,b 2 =0.01 












c=0.85,b 2 =0.002 














c=1,b 2 =0.002 . 




















c=0.85,b 2 =0.01 


i 
1 












c=0.85,b 2 =0.002 

c=1,b 2 =0.002 







































-0=0.85, b =0.01 
-c=0.85,b 2 =0.002 
c=1,b 2 =0.002 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



0.2 0.4 0.6 0.8 



FIG. 3: a) Evolution of flu and ujd with the scale factor a. To plot these curves we have fixed the best-fit value of Q m o = 0.27. 
The solid, dashed and dotted lines stand, respectively, for the pairs (b 2 = 0.01, c = 0.85), (b 2 = 0.002, c = 0.85), and (b 2 = 0.002, 
c = 1). b) Evolution of H/H 2 with the scale factor a. Note that, as ujd is becoming more and more negative and crosses the 
phantom divide line (Panel 3a), the function H increases from negative to positive values. As in the previous Panel, the value 
of the matter density parameter has been fixed at Q m o = 0.27 and the solid, dashed and dotted lines correspond to the above 
combinations of the parameters b 2 and c. 



data 0, [U 



K = ^n^ -— = 1.70 ±0.03, (19) 
Jo E{z) 

where the dimensionless function E(z) = H(z)/Hq and 
zi s = 1089 ± 1. In order to obtain the distance, we need 
to find out the evolution of FId(z) and H(z), so we need 
to solve Eqs. (fT!?)) and (fT3")l numerically. Since the deriva- 
tives in Eqs. (fl~2f and (fl~3f are with respect to x = In a, 
we need to rewrite them with respect to z. We find 



dH 

dz 




V-3_l 
2 2 2 7 



(20) 



and 
dfl 



D 



Q 



D 



dz 



1 + z 



(l_^)(l + M^)- 36 2 + ^ 



(21) 

By solving numerically the above equations, we then ob- 
tain the evolutions of Qd and H as a function of the 
redshift. 



B. Results 

In Figs. (1) and (2) we show the results of our statis- 
tical analyses. Figure (la) shows the c—b 2 plane for the 
joinf analysis involving the 192 ESSENCE SNe la data [3] 
and the other cosmological observables discussed above. 
For this analysis the best fit values are Sl m o = 0-271q'q3, 
c = 0.85lg;J|, and b 2 = 0.002t°;°J 2 (at 68.3% c.l.) with 



yZ min = 196.16. If we fix c = 1, we find n m0 = 0.26±g;g| 
and b 2 = 0.005^;^ (at 68.3% c.l.) with x 2 = 198.96. 
The plane O m0 — c is shown in Panel (lb). Figure (2a) 
shows the same parametric space c—b 2 when the 192 
ESSENCE data is replaced by the new 182 gold sample 
0. In this case, we find Q m0 = 0.29 ± 0.04, b 2 = 0j£ 01 , 
c = 0.88i£;£! (at 68.3% c.l.) with* 2 = 158.66. The plane 
f^mO — c for this latter combination of data is also shown 
in Panel (2b). We note that from both combinations the 
value of b 2 is very close to 0, which suggests a very weak 
coupling or a noninteracting HDE. Such a result is also 
in agreement with the limits found in Ref . [l3| . By fixing 
b 2 = 0, we also fit the HDE model without interaction 
to the observational data discussed above. These results 
are summarized in Table [I] For the sake of comparison, 
we also list the best-fit results for a flat ACDM model. 

Finally, by fixing the value of the matter density pa- 
rameter at f2 m o = 0.27 we show, in Fig. (3a), the evolu- 
tions of r^£) and ujjj with the scale factor a. The behavior 
of H/H is also shown in Fig. (3b). The curves displayed 
in these Panels are complementary in the sense that from 
them, we see that while u>e> crosses the cosmological con- 
stant barrier to the phantom region, H increases from 
negative to positive values. The distinctive future super- 
acceleration, which is an evidence of a phantom behavior, 
is apparent from Panel (3b). 



IV. CONCLUSIONS 

We have analyzed the observational viability of an in- 
teracting HDE model characterized by the parameters 
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TABLE I: The best-fit results for the HDE parameters 



Model 


Results 


Gold 


Gold±A±ft 


ESSENCE 


ESSENCE±A±ft 




x 2 


158.27 


158.66 


195.34 


196.16 


With 
Interaction 


b 2 

C 


n on+0. 29 
u -°^-0.13 

n +0.2 

U -0 
fl 89+0-48 


0.29 ± 0.04 

n+0.01 
u -o 


r\ 97+O.23 
u - z ' -0.15 

0.85t°;« 


27+ - 04 

U-^ 1 _o.03 

002+ 001 
u.uuz_ 002 

8^+ ' 18 
u - oo -O.02 




x 2 


158.27 


158.66 


195.75 


196.29 


b 2 = 


C 


0.31+H 7 
0.82lg : « 


0.29 ± 0.03 

c = o.88iS;« 


n 9 7 +0.03 
u ' z '-0.14 

c = 0.85t°;« 


n 97+O.O3 
u - z ' -0.02 

85+ 01 
U.80_g 02 


ACDM 


x 2 


158.49 


161.87 


195.34 


196.12 






0.34 ± 0.04 


0.29 ± 0.02 


0.27 ±0.03 


0.27 ±0.02 



c and b 2 and shown that in order to make it compati- 
ble with current observational data only the future event 
horizon can be used as the IR cutoff. By assuming a 
particular interacting term and fitting the model to the 
observational data, we have found that an HDE is ca- 
pable of explaining the current observations of SNe la, 
BAO and shift parameter from CMB data for the inter- 
vals of the parameters fl m a, c and b 2 displayed in Table I. 
From both analyses (involving the 192 SNe la ESSENCE 
data and the new 182 Gold sample) we have also found 
that a very weak coupling or even a noninteracting HDE 
model is favored (b 2 ~ 0). By using the fitting results for 
the parameters, we have shown that the expansion rate 
in the interacting HDE model can be super-accelerated, 
i.e. H > 0, which is the key physical consequence of the 
so-called phantom behavior [l£|. This phantom regime 
is obtained after a transition from lud > —1 to lud < — 1 



(phantom crossing behavior) , which is a feature well stud- 
ied in the context of some non-minimally coupled scalar 
fields Q. 
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